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The general form of the global conservation laws for iV-body systems in 
the first post-Newtonian approximation of general relativity is considered. 
Our approach applies to the motion of an isolated system of N arbitrarily 
composed and shaped, weakly self-gravitating, rotating, deformable bodies 
and uses a framework recently introduced by Damour, Soffel and Xu (DSX). 
We succeed in showing that seven of the first integrals of the system (total 
mass-energy, total dipole mass moment and total linear momentum) can be 
broken up into a sum of contributions which can be entirely expressed in terms 
of the basic quantities entering the DSX framework: namely, the relativistic 
individual multipole moments of the bodies, the relativistic tidal moments 
experienced by each body, and the positions and orientations with respect to 
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the global coordinate system of the local reference frames attached to each 

body. On the other hand, the total angular momentum of the system does 

not seem to be expressible in such a form due to the unavoidable presence of 

irreducible nonlinear gravitational effects. 
04.25.Nx, 95.10.Ce 
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I. INTRODUCTION 



Recently, Damour, Soffel and Xu |l]-|J (DSX) developed a new, exhaustive approach to 
the first post-Newtonian dynamics of a system of N extended bodies. This theory is based 
on the complementary use of some local coordinate systems (attached to each body) and of 
a global coordinate system used to describe the orbital motion of the N bodies. Detailed 
analyses of the laws of global translational motion and local rotational motion P] of the 
bodies have been given. 

In this paper, we address the question of the global conservation laws of an iV-body 
system, and their link to the quantities introduced in the DSX framework. We use the 
word conservation laws to mean the first integrals related to the total four momentum and 
angular momentum of the system, and to the center-of-mass theorem. The existence, on 
the first post-Newtonian level, of these ten integrals is guaranteed by the general form of 
the field equations |5|-[10| . Our main problem concerns the form of these conservation laws. 
More precisely, we investigate whether they can be entirely expressed in terms of the basic 
quantities introduced in the DSX post-Newtonian theory: namely, the set of the relativistic 
individual mass and spin multipole moments of the N bodies (M^ , Sfy and the set of 
the gravitoelectric and gravitomagnetic tidal moments experienced by each body (Gf , H^). 
[Here, as in Refs. |T]-Q whose notation we follow, A,B = 1, N labels the various bodies, 
and L = is a multi-spatial index] The former fully characterize the structure of the post- 
Newtonian gravitational field generated by an extended body in its local coordinate system, 
while the latter characterize the tidal action of the other bodies in this local coordinate 
system, including inertial contributions due to its acceleration and rotation. Obviously, the 
positions and velocities, with respect to the global coordinate system, of the origins of the 
local systems attached to each body, as well as their orientations, need to be considered, 
and will also enter the final expressions. 

Let us start by discussing the various methods used for deriving some explicit forms of 
the global conservation laws in general relativity. First, let us remark that the expressions 
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based on surface integrals at infinity 0, [f|, |H| are of no real use within the post-Newtonian 
context because one loses a power 1/c 2 in reading out a conserved quantity from the asymp- 
totic behaviour of metric. Concerning approaches where the conserved quantities are given 
as three-dimensional integrals, we note the formulation of Fock || , taken up by Brumberg 
TT| , and that of Chandrasekhar and coworkers |7J , || , followed by the parametrized version 



of Will [ 12| . Neither approaches constitute useful starting points for us. Indeed, on the one 



hand, Fock and Brumberg restricted themselves to the (ill-defined) case of a "rigidly rotat- 
ing" bodies and introduced some mass moments which are not compatible with the ones 
that enter naturally the DSX formalism, while, on the other hand, Chandrasekhar et al. 
restricted themselves to the special case of perfect fluids and chose basic variables which do 
not fit well within the DSX scheme. 

A more convenient starting point for our purpose is the work of Blanchet, Damour and 
Iyer [J3], [|TJ|] (see also |15|)- These authors have defined global-frame post-Newtonian mass 



and spin moments of an arbitrary, isolated system of bodies in the form of integrals over the 
compact supports of the bodies, and, they explicitly checked that the lowest moments were 
conserved quantities. Our task here will consist in transforming their expressions, involving 
integrals over the global-time simultaneity surface, into a combination of terms involving 
integrals over N separate local-time simultaneity surfaces. Notice that nothing guarantees 
a priori that these manipulations will lead to final expressions containing only the "good" 
moments introduced in the DSX scheme. In fact, we shall succeed in this task for seven of 
the globally conserved quantities and fail for three of them (the three components of the 
trickier total angular momentum). 

Let us note in advance that there are several limiting cases for which it is already known 
that some of the globally conserved quantities can be entirely expressed in terms of some dy- 
namically relevant individual multipole moments. First, there are the cases where one trun- 
cates the multipolar series: keeping only the monopole contributions defines the "Lorentz- 
Droste-Einstein-Infeld-Hoffmann" model (LD-EIH), while keeping also the intrinsic spins of 
the bodies define the "pole-dipole" model (PD). In both cases, we dispose of well established 



forms of the conservation laws (see e.g. ||, for the monopole case, and [I(J, [17 for the 
pole-dipole case). These will provide useful checks on our general results. Another limiting 
situation of possible relevance is that of a test point-like mass moving in the gravitational 
field of iV — 1 extended bodies. The motion of the artificial/natural solar system satellites 
is a typical case of this category. Ref. [Q computed explicitly the form of the corresponding 
Lagrangian in terms of the gravitational potentials (W, W a ), which can be algorithmically 
constructed using the formulas given in the Appendix of ||. However, in this case we have 
first integrals only when the Lagrangian (describing geodesic motion) possesses some contin- 
uous symmetries (Noether's theorem). For instance, we can consider N = 2 with Mi << Mi 
(restricted two body problem) and with a stationary and/or axially symmetric central body 
(see e.g. Ref. [0 for a study of a restricted problem of this type). This type of limiting 
situation will not give us useful checks. 

Finally, let us remark that we hope that the present work will find practical applications 
in the relativistic motion of binary stars or in the celestial mechanics of the solar system. 
Let us recall, for instance, that one simplifies the dynamical ephemeris of the solar system 
(at the Newtonian approximation) by eliminating the motion of the Sun via the algebraic 
relation expressing that the total dipole mass moment of the solar system vanishes for all 
time (in a suitably mass-centered frame; [|P9| ). 

In Sec. II we consecutively discuss the cases of the mass-energy integral, linear momen- 
tum integral and briefly comment on the angular momentum integral. Sec. Ill contains a 
summary of our results. Throughout the paper we follow the terminology and notation used 
in Refs. [HI 



II. FIRST INTEGRALS OF THE DYNAMICAL LAWS 
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A. Comments on the change of time-simultaneity integration domains 



As mentioned in the preceding section, transformation of integrals from the global to 
several local time simultaneity surfaces is a common point to all particular cases of conser- 
vation laws. We shall thus start our discussion with a brief technical comment concerning 
such transformations. 

Considering a specific body A, we study some integral 

l A (t) = ( d 3 xf(t,x), (2.1) 

J A 

performed on the global time t = const, surface, where f{x^) is some given function defined 
on the compact support of body A. We seek a transformation of the right hand side of X A 
such that it can be written as an integral on a local-time T A = T A (t) = const, simultaneity 
surface, say 

T A (T») = / d 3 X A F(T° A (t),X A ) . (2.2) 

J A 

Here, T A (t) denotes the value of the local time (in the reference system attached to body A) 
corresponding to the event on the central worldline of the body-v4 reference system (X^ = 0, 
also quoted as C A ), whose global time coordinate is equal to t. In equations, if we write 
the spacetime coordinate transformation between the global coordinate system x M = (ct, x l ) 
and the body- A local one X A = (cT A , X A ) as 

x» = a?(X% = z» A (T A ) + e» Aa (T A )\X a A 

+^A A X\ - ±(A A .X A )X a A ] , (2.3) 

T A {t) is defined as the unique solution of ct = x°(T A , 0), i.e. ct = z A {T A {t)). 

First, let us denote by f A (T A ,"K A ) the original function f(t, x) reexpressed in terms of 
the local spacetime variables X A : f A (X A ) = f{x^{X A )). By mathematically transforming 
the variables of integration in Eq. Q2.1|) we get 

(3) \ - 1 

f A [T A (t,X A ),X A ] , (2.4) 

t=const J 
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where \dX/dx\f2 CO nst is the spatial Jacobian det(dX a / dx l ) (a,i = 1,2,3) computed when 
keeping fixed the value of t, and where T A (t,X A ) denotes the solution of ct = x°(T A , X a ). 
From Eq. ( |2.3|) (or Eqs. (A5) of Ref. 0]), the latter quantity reads explicitly 

T A (t,X A ) = T° A (t) - \v%X\ + 0(4) , 

so that by expanding f A [T A (t, X^), X^] in powers of the small time shift (V^.X^) /c 2 = 0(2) 
we can express it within a sufficient accuracy in terms of functions computed on a local-time 
simultaneity surface, namely T A = const. = T A {t): 

f a [T A {t, X A ) , X A ) = f A (T° A (t), X A ) 

-±VZX a A d T f A (T° A (t),X A ) 

+0(4) . (2.5) 

As for the three-dimensional Jacobian entering Eq. (|2.4j ) it is easy to see that it can be 
expressed as 



dX 


(3) 


/ dt N 




dX 


(4) 


dx 






X a =const 


dx 





e° (T°) + -A A X a 



dX 



dx 



(i) 



+ 0(4) . 



(2.6) 



where \dX/dx\^ = det(dX a / <9x M ) is the full four- dimensional Jacobian associated with the 
coordinate transformation ( |2.3| ). The time derivative (dt/dTA) in Eq. ( |2.6| ) is obtained from 
Eq. ( |2.3|) or from putting Vs = in expression (A6) in Appendix A of [[|. 
Finally, we get 



rf 3 A/ 



<9X 



<9x 



(i) 



-i 



e % (^a) + ~^-^A^A 



/a (T a , Xa) -V A X A d T fA (T A , X a 



(2.7) 



where it is convenient to leave unexplicated the four- dimensional Jacobian because it will be 
directly cancelled when using the transformation laws of the mass and mass current densities 
<7 M entering f A (t,x). For completeness, let us however mention its value 

7 



= l-^-W"(T,X) + 0{4) , (2.8) 
& 

where W"(T,X) = G" A - A\X a + 0(2) is the inertial contribution to the local potential 
due to the change of the time scale and the acceleration of the body A frame. Here, G" A = 
c 2 \n(dT I drf) a = v\/2 — c 2 \ne A0 + 0(2) measures the relative scaling, along the central 
worldline X A = 0, between the local time T A and the global Minkowskian proper time 
dr f = J-f^dz^dz'/c (see Sec. VI.E of [§). 

Note that, in geometrical terms, the mathematical transformations we have just per- 
formed correspond to using a mapping between the t = const, and T A = const. = T A {t) 
hypersurfaces by means of the congruence of worldlines Cx^ of constant spatial local coor- 
dinates (see Sec. III.D in Jl|). 



dX 
dx 



B. Mass-energy integral 



The Blanchet-Damour post-Newtonian total mass-energy m(t) of an isolated system can 



be written as [13 



m(t) 



d 3 x a(x) -~ I d 3 x (aV) + 0(4) 
c z dt 



(2.9) 



where <7 M are densities of mass and mass currents [a = (T 00 + T u )/c 2 and a 1 = T 0l /c; T^ v 
denoting the components of the stress-energy tensor in the global coordinate system]. The 
integration in (|2.9|) is to be performed over a global-time t = const, hypersurface spanning 
the whole iV-body system. However, as is nonzero only in the neighbourhood of the 
bodies, we can directly use the results of the previous sub-section to decompose m(t) as a 
sum of iV terms integrated over local simultaneity surfaces T A (t) = const.. 

Let us recall the transformation law between the global and local coordinate systems of 
the mass densities pertaining to a given point of body A 

OX 



a(x\ 



dx 



(4) 



l + 2^fW) 



+-VXZ a A (X) +0(4) 



(2.10) 



as given in |L]. 

Putting together the definition ( |2.9| ) and the method explained in the previous sub- 
section II. A [where we note that the four dimensional Jacobian cancels between Eqs. (p.7[) 



and (|2.10|)1 we arrive after some algebra to 



m(t) = £ 

A 



. \ ( i 1 . .2\ , 1 ^ / i fir •« 



M [ 1 + ^) + ^Ht a ( M ' V ' 



+ 0(4) , 

(2.11) 



where all quantities on the right-hand side must be evaluated at the intersection of the 
t = const, hypersurface with the central worldline of the corresponding body (i.e. for 
T A = T%(t)). We recall that the quantities M A (which are not constant in general) denote the 
local, individual relativistic mass monopoles of each body. Each M A is a directly observable 
quantity in the sense that it is just the gravitational mass measured from interpreting the 
locally measured orbital motion of artificial/natural satellites around body A. Similarly 
the M A (for I > 1) are the locally measured mass multipole moments of body A, and G A 
the locally felt tidal moments. The other quantities entering Eq. ( j2.11 ) are related to the 



way the local A-reference system is moving with respect to the global coordinate system [in 
particular G'[ A = (G" A , — A A , +3A A a A A > /c 2 , 0, 0, ...) measure the inertial contributions to 
the tidal moments felt in the local A system] . The result ( p.ll|) is new, and its precise form 
(e.g. the numerical factor 21 + 1 in front of the M A G A product) is different from what one 
might have naively expected from the standard Newtonian expression for the total energy 
(e.g. fU). 

As discussed in detail in Ref. [J] the DSX framework leaves open some freedom in fixing 
several quantities related to the origin an orientation of the local coordinate systems, as 
well as the gauge for the time coordinate along the worldline La- We shall call "standard 
worldline data" the case where this freedom is used to satisfy the following constraints: (i) 
WA , VTa , M a (Ta) = (which means identifying the origin of all local coordinate systems 
with the relativistic mass centers of the bodies), (ii) VA , VTa , W a (Ta, 0) = (the so called 
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weak effacement condition for the external gravitational potential in the local frames). Note 
that we still leave unconstrained the orientation of the local frame axes which can undergo 
a slow rotation described by the matrix R l Aa (T). In the case of the standard worldline data, 
Eq. (p.lip simplifies to the form 

1 



'(*) = £{m 

A l 



standard f+\ \ * J j\ j A 



+ ^7^^ + 0(4), (2.12) 



2c 2 ^ l\ 

where G' A = Eb^aG b/a = Eb^aW B (z a ) + 0(2) [G' A = -G" A for standard data] denote 
the value on the central worldline La of the Newtonian potential generated by all the other 
bodies [w B (z A ) = Ei>o ( -^GM*d£\z A - zb^ 1 }. 

In the monopole (LD-EIH) or pole-dipole limit Eq. (|2.12|) yields the well-known result 
that the total mass-energy is the sum of the total rest-mass and the kinetic and potential 
(oc GM A M B /\za — zb\) energy terms. Let us recall that, in the general case, the individual 
post-Newtonian gravitational masses M A are no longer constant (because of tidal forces 
acting on the extended bodies, see Eqs. (4.20a) and (4.21a) of 0) and that it is not a priori 
evident that the quantities m(t) or m standard (£) are conserved. As a check on our algebra, 
we have verified by a direct calculation that this is indeed the case. 

As an aside, let us conclude this sub-section by noting that if one defines, as an auxiliary 
technical quantity, the "Fock mass" of the A-th. body by the relation c 2 Mp = J A d 3 X A (1 + 
W/2c 2 )T A ° (where denote the X A -X° A component of the stress -energy tensor considered 
in the local coordinate system X A ), one can write our result ( ]2.12j ) in a formally compact 
(and familiar looking) form. The name we give to this quantity is based on the fact that 
in his book |J Fock used such an expression for the total mass-energy. Note, however, 
that he always used it in the global coordinate system. Our definition is written in the 
local system X A , but we use for W(T, X) the total potential in the local frame, containing 
both internally and externally (including inertially) generated contributions. By using the 
expressions (4.15) of 0] for the tidal expansion of W (T, X) we find (in any worldline gauge) 
the following relation between the Blanchet-Damour mass and the Fock one: 
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M 



F 2c 2 4" /! 



M^Gl + 0(4) 



(2.13) 



Thus in the case of standard data we can rewrite our previous formula (|2.12|) in the following 
form 



m 



, standard 



(*) =J2 M F 

A 



1 + 



2c 2 



(4 - G'a) 



0(4). 



(2.14) 



Let us, however, emphasize again that it is only the Blanchet-Damour mass M A which is 
directly observable as a gravitational mass determining the orbital motion of satellites of 
body A. The Fock mass Mp is just a mathematical construct. 



C. Center-of-mass integral and linear momentum 

The Blanchet-Damour post-Newtonian dipole mass moment rrii(t) of the whole system 
13(1 can be written as 



rriiit) = / d 3 xx l a(x) 



(2.15) 



-~ J d 3 xai (x*x j - U l3 x 2 ^ + 0(4) . 



It satisfies the following conservation law |13 

d 2 rrii(t) 
dt 2 : 







(2-16) 



In fact, the first derivative of rrii(t) is nothing but the conserved total linear momentum of 
the system 

drrii(t) 



Pi 



dt 



const. 



(2.17) 



Employing the results of Sec. II. A and definition ( |2.15| ) we obtain after tedious but 
straightforward calculations the following form of the total mass dipole moment rrii(t) of the 
system 
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rrii 



1 ... "X 
^ V A S A ~ ^ a A m A 



where s l A 



+ X/ vi 

A 

+0(4) 



e ijk s A 



l + ^(~v 2 A + G A 



Air 

a " a 



(2.18) 



?-Lt A„>J A, Hi A 



e Aai e A 



a2 ...eX n M a / a2 - an , where M^ 1 ) 



dM A /dT A and where the brackets in R A [ a R Ab ] mean antisymmetrization [it^i^] = \{u a Vb 



Uf,v a )]. It should be emphasized that the fact that the final expression ( 2.18 ) can be entirely 
written in terms of the "good" moments entering the DSX framework is far from being 
a trivial result. In the intermediary calculations the "bad" moments (Nl, Pl) defined in 
Eqs. (4.22) of |2| enter at several places before finally cancelling. We also remark that the 
last sum in curly brackets vanishes if one uses standard worldline data, as defined above. 

As mentioned in Sec. I, we can get partial checks on our results by considering models 
where the multipole series is highly truncated. In particular, if we keep only the mass 
monopoles of the bodies (LD-EIH limit), formula ( |2.18| ) reduces to 



m 



LD-EIH 



A 



A J, 



1 

2c~ 2 



l + ^(v\-G' A 



A J. 



1 



B^A 



+ 0(4) 

GM B 
tab 



+0(4) , 



(2.19) 



where the second row applies to the case of standard worldline data. Eq. (|2.19| ) agrees with 
previous results 0, flOfl . In the case of the pole-dipole truncated model the expression for 
the mass dipole reads 

1 



m 



PD 



m 



LD-EIH 



(2.20) 



a result previously derived by Damour and Schafer |J7] from the spin dependent Lagrangian 
of Ref. |T§. 
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Let us remark, as an aside, that defining some "Fock" local mass dipole moments for 
instance by c 2 M A F = j A d 3 X A X A (l + W/2c 2 )T 00 does not simplify at all the writing of our 
result (|2.18| ). In fact, this definition introduces several bad algebraic structures (notably 
the moments N A ; see which do not enter the final dynamical results of the DSX 

formalism. This is one of the reasons why Fock, in his book |5]], did not succeed in getting 
a good definition of the mass centers of the individual bodies (in spite of the fact that, in 
the case of the mass center of the entire isolated system, the Blanchet-Damour and Fock 
definitions, written in the global coordinate system, give the same result; see Eq. (3.45) of 

m 



From the result ( |2.17| ) above, we can easily derive the following explicitly DSX-like 
expression for the total linear momentum of the iV-body system: 



1 + 



c 2 \2" A 



v\ + G" A 



1 



3 M 
~2 U 'A° A 



3 d 



-~A E \ [lMi {l) G A L > +(/ + !) M£G A L W) - i {v\z\ - 5^ A .z A ) R Aa £ ±M*G; 



iAl 



~ 72 ( Z A Z A ~ 7^i Z l) RAa E T\~rFT\ {^lG A l) 2 Z A^A[ a RAb] E 



j" l\ dT A V L aL l c 2 A[a Abi Y 11 dT A 

K) (i + + >4 «<) + ^ A ±f A - 1 

+0(4) . 



v A m A a A 



X21) 



Again, the last sum of terms in curly brackets disappears in the standard worldline gauge. 

The LD-EIH form of the linear momentum is obtained by retaining only the mass 
monopole terms 



Pld-eih 



eK 



2c 2 r B % rAB 



G v GM A M B 



P-AB-V B )n AB 



-0(4) . 



(2.22) 



The presence of spin dipoles (in the framework of the PD model) results in the following 
additional term 
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PpD = P'lD-EIH + 4 E «M( + 0(4) , (2.23) 



r 2 



as mentioned in IT? 



As in the case of the total mass-energy, but with more work, one can directly check that 
the global-time derivative of Eq. (|2.21|) vanishes to the indicated accuracy. 

Due to the incompatibilities between the DSX and the Fock approaches mentioned above, 
it is not possible to compare directly our results with those given by Fock || or Brumberg 

HI. 



D. Angular momentum 



Damour and Iyer JL4| and Ref. [|3] have shown that the treatment of the local individual 
spin of bodies, members of an interacting A-body system, faces serious problems due to 
the unavoidable intervention of nonlinear, interbody gravitational effects. Mathematically, 
this complication manifests itself through the occurrence of "bad" DSX moments (Nl,Pl)- 
It has been, however, possible to reach a successful formulation of the individual rotational 
laws of motion through a carefully adjusted definition of the individual spin of each body |§. 
In the following, we briefly address the problem of breaking up the total angular momentum 
into a sum of contributions which make sense within the DSX framework. 

1, 0, I) 



The global angular momentum of the system reads 



Si{t) = e ijk / d 3 xx 3 \ a 



a 

72 



4 ' 
1 + -rw 



4w fc + -d k d t z(x, t) 



+ 0(4) 



(2.24) 



where z(x,t) = G J d 3 x' o[x! , t)|x — x'|. For an isolated system of arbitrary bodies it has 
been previously shown that Si(t) is conserved at the post-Newtonian level. 

In order to preserve the post-Newtonian accuracy [i.e. modulo 0(4)} of the result, one 
needs the transformation law of the global mass current a k (x) to local coordinate quantities 
with corresponding precision. After some algebra one arrives at 
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a [x] 



dX 


(3) ^. 


dx 





V A 6 /InS 



1 r 



+ - 



+- 



^ (**V*E) + P^ T (X 6 VJ5? 



-\ab 



S ab T cc 



+0(4) , 



^2.25) 



(where = de Aa /dT A , etc.) generalizing the formula cr J = t>^£ + P^ a £ a + 0(2) used 
throughout the series of papers [Note the three-dimensional Jacobian, as defined in 

Eq. above, in front of Eq. (ggg)]. 

Inserting this relation into the defining integral Q2.24j) and using the method of Sec. II. A 
one obtains an expression of the form 



Si (*) = £ [*f (M L , S L , G L , H L ; P L , N L ) + 9i + 6, 



(2.26) 



where tyf is a function of both the good DSX moments (Mr,, Si, Gl, Pi) and the bad ones 
(Pi, iVi), and where the remaining two terms (which follow from the last term in Eq. ( |2.25| )) 
read 



e ljk z\VXR k M J A d 3 X A 

° 2 ^2i — e ijkR\oYA Rk Ac I A ^ X A X 



ia6 



S ab T cc 



T bc - 5 bc T dd 



Note that if we augment our list of "bad" moments by including Q^ b = J A d 3 X A T ab and 



Q 



ab:c 



J A d 3 X A T ab X A , we can express the total spin Sj(t) in terms of some individual 



moments of the iV-bodies. We tried to get rid of the non-dynamical moments (P^, N£, Qfy 
by using the local conservation of energy-momentum (V a T at3 = 0) to connect them to the 
dynamical moments (M^,S^) and the tidal ones (G^,H^). We did not succeed in doing 
so. [In fact, such transformations introduced an undesirable dependence upon the internal 
part W2 + (X) of the local gravitational potential]. We thus hypothesize that the total spin 
Si of the system is not algebraically reducible to the DSX dynamical quantities. 
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III. CONCLUSION 



The algebraic form of the global relativistic conservation laws has been examined within 
the perspective of the DSX post-Newtonian dynamics of an A-body system. We succeeded 
in breaking up seven of these conservation laws (mass-energy, center-of-mass quantity, and 
linear momentum) into a sum of individual contributions involving only the basic dynam- 
ical quantities of the DSX formalism. The angular momentum conservation law resisted, 
however, our efforts. 
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